Electron Shielding 


In the above laws, there is one additional factor which must be in- 
serted. Surrounding all the ions in the center of a star is a cloud 
of electrons (and some positrons); the net negative charge of this 
cloud alters the Coulomb potential each ion sees and thus increases 
the reaction rate. To estimate the importance of each term, first 
consider the Boltzmann equation, which describes the relative den- 
sity of particles. If the statistical weights of two states are the same, 
then 
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where F is the energy difference between the two states. Now sup- 
pose there is a varying electrostatic potential, (r). For normal 
stars, Zed < kT, so we can write the density of particles, as a 
function of @ as 
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If we separate out the ions and electrons, then 
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and the charge density as a function of the potential is 
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This can be simplified by noting that the average charge density 
must be zero, t.e., 
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When we substitute in the definition for n; and ne via (5.1.2) and 
(5.1.4), this becomes 
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where x; is the mass fraction and f; is the ionization fraction. Since 
the conditions are such that nuclear reactions are occurring, f; = 1, 
so the charge rey finally becomes 
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This charge distribution now allows us to solve for the potential 
explicitly, (at least in the spherically symmetric case) via the Poisson 
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l d do 
2 2 
— LG | 
wg r2 dr (r N oe 
The solution to this equation (which you can check if you like) is 
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where rp is the Debye-Hückel length 
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From the above equation, the potential goes to the normal electro- 
static potential, i.e., 6 — Ze/r as r — 0, so the electron shielding 
only works effectively at large distances, i.e., r > rp. However, most 
reactions are due to particles near the Gamow peak; the classical 
separations of these particles are rg = ZaZxe?/Eo. A comparison 


shows that 
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For all normal stars, rp >> rg, which means that electron shielding 
is not very effective. Electron shielding changes the potential barrier 
(11.3) slightly, by 
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and thus increases the cross sections and reaction rates by f = 
e£D/kT with 
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For typical stellar densities and temperatures, this correction will be 
~ 10%. This “weak screening” approximation, however, does not 
apply to the high density-low temperature regime of white dwarfs. 
Under those conditions, the form of the equation will be different. 


Summary for Reactions 


When electron shielding is included, the nuclear reaction rate equa- 
tion becomes 
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or, in terms of number density 
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For non-resonant reactions, (ov) can be computed directly from 
S(Eo), Caen Za, Zx, Aq, Ax, and T, while for resonant reac- 
tions, (ov) is derived from (wy), Ep, and T. 


Now let’s simplify the terminology by first assigning 
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and then defining the lifetime of species X against reactions with 
particle a to be 
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The form of this equation is now identical with that associated with 
radioactivity or other decay phenomena. (Unlike radioactive de- 
cay, of course, Ta depends on the external environment.) From the 
equation above, 
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Note that the Kronecker delta disappears from the final lifetime 
calculation. According to (12.2.1), the reaction rate for identical 


particles has a factor of two in the denominator, but since each 
reaction destroys two particles, this factor cancels out. Note also, 
that by (12.2.4) and (12.2.5), the total lifetime of particle X from 


all reactions is j i 
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